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Rescue (SAR) System (COINSS). This system provides several important features not 
provided by current inland SAR computer systems. First is the ability to model 
movement of the target. Second is modeling the effect terrain has on the movement of 
the target. Third is the visual presentation of a probability map, a color display showing 
the probability that the target is located at various geographic positions. COINSS is 
developed in the Java programming language. It is designed to be implemented with a 
map-based planning system using loosely coupled components. COINSS provides the 
initialization, movement, and search algorithms which are used by the planning system to 
support the search operation. The initialization algorithms define the search area where 
the SAR operation will occur. Initial areas are defined for the target. COINSS models 
the movement of the target as a discrete time Markov chain. Bayes theorem is used to 
update the probability map when negative search information is provided. This thesis 
will improve inland SAR operations by providing the first model with an interactive 
graphical user interface and a model of target movement. 



v 




- 



VI 



TABLE OF CONTENTS 



I. INTRODUCTION 1 

A. BACKGROUND 1 

1 . The Responsibility of Land-Based Search and Rescue 1 

2. SAR Tools and Process 1 

B. PROBLEM 3 

C. APPROACH 4 

1 . The Key Elements 4 

2. Two Concepts of Movement 6 

D. COINSS 7 

E. THESIS STRUCTURE 8 

II. THE SEARCH AND RESCUE PROCESS 11 

A. INITIAL ACTION 11 

1. Obtaining an Initial Position 1 1 

2. Defining the Search Region 12 

B. SEARCH ACTION 12 

1. Identifying the Search Area 13 

2. Search Execution and Results 13 

III. LOGICAL MODEL 15 

A. PROBABILITY MAP COMPONENTS 15 

1. Elements of the Grid Square 15 

B. INITIALIZING THE PROBABILITY MAP 1 6 

C. UPDATING THE PROBABILITY MAP 17 

1. Time-Based Update 17 

2. Search-Based Update 23 

IV. VISUAL REPRESENTATION OF THE MODEL 25 

A. COMPONENT-BASED DESIGN 25 

B . BASIC CONSTRUCTION OF THE VIEW 25 

1. Grid Dimensions 25 

2. Probability Display 26 

C. INTERACTING WITH THE VIEW 27 

1. Selecting a Geographic Map 27 

2. Defining the Search Region 28 

3. Defining the Initial Areas 28 

4. Defining the Search Area 29 

V. CONCLUSIONS AND RECOMMENDATIONS 31 

APPENDIX A. MARKOV MOVEMENT MATRICES 33 

APPENDIX B. COINSS SCENARIO 39 



vii 



APPENDIX C. COINSS CODE 55 

LIST OF REFERENCES 81 

INITIAL DISTRIBUTION LIST 83 



Vlll 



EXECUTIVE SUMMARY 



The National Search and Rescue Plan establishes the United States Air Force as 
the executive agent for Inland Area Search and Rescue (SAR). This area encompasses 
the Continental United States, except Alaska, and waters under the jurisdiction of the 
United States. The Air Force established the Air Force Rescue Coordination Center 
(AFRCC) as the contact point for federal resources for inland SAR [18]. 

In addition, the National Search and Rescue Plan dictates that Search and Rescue 
Coordinators should provide for the fullest cooperation and use of additional Federal 
military, civil, State, local, and private agencies as may be necessary and practicable [18]. 
The On-Scene Commander is responsible for overseeing the execution of the entire 
search plan. The Commander uses search tables to formulate each assignment in the 
overall search plan. The individual search assignments are then provided to each search 
unit for implementation. The results of the individual search assignments are 
communicated back to the Commander and used to update the search scenario. Based 
upon the updated scenario, the Commander reevaluates the search plan, creates new 
search assignments, and the process is repeated. 

Currently, numerous diverse software tools and spreadsheets are available to help 
do numeric planning calculations and table referencing. These stand-alone aids compute 
the data for the Commander, who must then interpret the data and incorporate it into the 
search plan. Additionally, management programs are available for ongoing analysis of the 
search plan. 
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However, despite the availability of a variety of search software programs, three 
key factors prevent their widespread use and limit their usefulness during an active inland 
search operation. These factors are: 

1. the lack of a graphical mapping interface, 

2. the lack of inter-operability among the computer systems of cooperating 
agencies, and 

3. the inability to account for a mobile target. 

This thesis designs and implements the Coordinated Inland Area Search and 
Rescue (SAR) System (COINSS). This system accounts for the three factors not 
provided by current inland SAR computer systems. First, COINSS provides a color 
display of the target probability distribution, which can be used as an overlay on a map- 
based planning system. Second, it is developed in the Java programming language and is 
designed to be implemented with a map-based planning system using loosely coupled 
components. Finally, it models the movement of the target as a discrete time Markov 
chain. 

COINSS provides the initialization and search algorithms which are used by the 
planning system to support the search operation. The initialization algorithms define the 
search area where the SAR operation will occur. Initial areas are defined for the target. 
Bayes theorem is used to update the probability map when negative search information is 
provided. 

This model will improve inland SAR operations by providing the first model with 
an interactive graphical user interface and a model of target movement. 
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THESIS DISCLAIMER 



The reader is cautioned that the computer programs developed in this research 
may not have been exercised for all cases of interest. While every effort has been made, 
within the time available, to ensure that the programs are free of computational and logic 
errors, they cannot be considered validated. Any application of these programs without 
additional verification is at the risk of the user. 
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I. INTRODUCTION 



A. BACKGROUND 

1. The Responsibility of Land-Based Search and Rescue (SAR) 

The National Search and Rescue Plan establishes the United States Air Force as 
the executive agent for Inland Area Search and Rescue (SAR). This area encompasses 
the Continental United States, except Alaska, and waters under the jurisdiction of the 
United States. The Air Force established the Air Force Rescue Coordination Center 
(AFRCC) as the contact point for federal resources for inland SAR [18]. 

In addition, the National Search and Rescue Plan dictates that Search and Rescue 
(SAR) Coordinators should provide for the fullest practicable cooperation and use of 
additional Federal military, civil, State, local, and private agencies as may be necessary 
and practicable [18]. To this end, the AFRCC has established agreements with each 
State, which describe the responsibilities for different emergencies. However, these 
agreements are not consistent from State-to-State, and often do not fully address the 
availability of local and private agencies to assist with the search effort. 

2. SAR Tools and Process 

Search planning is a well-researched and documented discipline. The National 
Search and Rescue Plan provides a complete set of tables and charts for creating and 
calculating search areas and search patterns. 

The On Scene Commander is responsible for overseeing the execution of the 
entire search plan. The Commander uses search tables to formulate each different search 
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assignment in the overall search plan. The individual search assignments are then 
provided to each search unit for implementation. The results of the individual search 
assignments are communicated back to the Commander and used to update the search 
scenario. Based upon the updated scenario, the Commander reevaluates the search plan, 
creates new search assignments, and the process is repeated. 

The bulk of search planning is still done on paper charts and maps. The 
Commander must take into account a large range of factors such as the target's last 
known position, time of that position, the terrain that is being searched, the type of target 
that is being sought, the number, experience, and composition of the available search 
units, etc. The Commander uses this information, as well as personal experience and all 
available charts and tables, to lay out the search area, divide the area into appropriate 
search sectors, assign individual search teams to cover the sectors, and establish the 
search pattem(s) that will be used by the search teams. 

Currently, numerous diverse software tools and spreadsheets are available to help 
do numeric planning calculations and table referencing. These stand-alone aids compute 
the data for the Commander, who must then interpret the data and incorporate it into the 
search plan. 

Additionally, management programs, such as CASIE III [3], are available for 
ongoing analysis of the search plan. These programs allow the searcher to define a 
search region and divide that region into search areas. Each area is assigned a probability 
that the target is located within that area. The searches are then conducted in the various 
search areas and the results used to update the probabilities by applying Bayes theorem. 
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An unsuccessful search in an area reduces the probability that the target is in that area. 
This reduction in probability is directly related to probability of success of the search 
method used. 

B. PROBLEM 

Despite the availability of a variety of search software programs, three key factors 
prevent their widespread use and limit their usefulness during an active inland search 
operation. These factors are: 

1 . the lack of a graphical mapping interface, 

2. the lack of inter-operability among the computer systems of cooperating 
agencies, and 

3. the inability to account for a mobile target. 

CASE III is an example of current software that provides the Commander with 
the raw data required to make search decisions. The program takes user input search 
areas and target location probabilities and the program provides, as output, probabilities 
that the target is located within the predefined search areas. When the results of a search 
are entered, the program then updates the probabilities, and outputs these revised values. 
The Commander is then responsible for translating that information to a physical map 
representation. Then the Commander assesses the search region and evaluates the next 
step to be taken in the search execution. The new search assignments, and eventually 
search results, must then be entered back into the program, where the cycle resumes. 
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The second factor that must be considered is the wide assortment of resources that 



can be brought to help execute a SAR case. The range of available personnel resources 
goes from trained, full-time federal and state professionals all the way to sparsely trained 
local volunteer organizations. The computer resources that will be available from these 
sources are unpredictable at best. In order for a system to be useful for all organizations 
participating in the search, the program must be able to run or be viewed on any 
computer system. 

The third factor is the abilities of the program itself. No model currently exists 
which accounts for a moving target on land. Probabilities are assigned to the various 
areas of the search region. After these probabilities are assigned, the target is assumed to 
be immobile. The changes in the probabilities in each area are only the result of 
unsuccessful searches that were conducted. 

C. APPROACH 

This thesis develops the Coordinated Inland Area SAR System (COENSS) to 
provide one method of resolving the problems addressed above. The design of this 
computer program was based upon several key elements. 

1. The Key Elements 

a. Portable 

Due to the urgent and unpredictable nature of SAR, it is not possible to 
prepare in advance for all possible scenarios. The Commander must make best 
use of the resources available at the time. The computer program must be small 
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enough to be easily and quickly loaded from a disk or downloaded over a network 
and provide a user interface that requires little user training. Additionally, the 
program should be able to operate on a wide variety of systems. 

b. Graphical interface 

A graphical interface would allow the Commander to do the SAR 
planning on the system, and not require a transfer of information between a paper 
chart and the computer model. The visual display of probabilities would provide 
a better description of the search progress. Additionally, the direct interface 
between the visual display and the search logic would provide a quicker update to 
the search execution once other search information has been received. 

c. Interoperable with available maps 

Again, SAR is unpredictable. We cannot know where the next SAR case 
will occur. The program should allow for use of maps from any electronic source 
to provide the graphical planning interface. 

d. Allow for target movement 

A target may not remain stationary. Ignoring the motion of the target in 
the execution of the SAR plan significantly impacts the location of the search 
areas, the search method used in those areas, and the results of the search. 
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e. Common picture view 



One advantage of the portability element is that all units coordinating on a 
search would be able to view the planning tool to assist in their search execution. 
The Commander would maintain and update the operational view of the system, 
which would be distributed to each search unit. 

2. Two Concepts of Movement 

Two choices are available to describe the target movement through the search 
area. Wagner [19] explains one method as being represented by a bundle of simulated 
tracks that encompass the many possibilities for target initial position and movement. 

The number of necessary tracks depends upon the complexity of the target being 
modeled, and should be set to a sufficiently large number in order to cover the complete 
range of possible initial conditions. A probability is associated with each track. This is 
the probability that the target is following the respective track. Wagner further explains 
how the tracks are monitored and updated, based upon negative search results. These 
tracks move over a grid search area, and require the use of random number generation to 
determine the progressing course of the simulated tracks. 

The second method is a pure stochastic representation of movement [19]. A grid 
is placed over the search area. Each grid square, S x>y , contains a velocity array, V^y,^, 
which contain elements that are the probability that the target is located in the square, 
moving with velocity i in the east- west direction, and velocity j in the north-south 
direction. The total probability for each square is computed by: 
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During time-based updates, the probabilities are moved from square to square, 
based upon their velocity, and then the velocity state is updated based upon a Markov 
chain. Bayes theorem is used to update the probability distributions when negative 
search information is received. This method requires no random number generation. 

D. COINSS 

COINSS is designed as a component to be used in conjunction with a map-based 
planning system. The map-based planning system must have certain capabilities and 
provide distinct information to COINSS. First, the map-based planning system must be 
able to provide positions on the map by latitude and longitude. This is required when 
initially creating the search region, as well as when initial areas and search areas are 
created. Second, it must also be able to provide the number of pixels, both horizontally 
and vertically, between two points on the map. This is required when initially creating 
the search region. Finally, it must be able to receive numerical inputs that represent 
probabilities. This is required for initial area probabilities and search area probability. 

COINSS is written in the Java programming language and is platform 
independent. This means that compiled code can be executed on a variety of computer 
systems that have the Java Virtual Machine (JVM), for example, Windows 95, 98, NT, 
and various UNIX platforms. COINSS can be loaded to any computer with the JVM and 
executed immediately. The map-based planning system used to demonstrate COINSS is 
also written in Java. Therefore, the total system can be distributed to execute on any 
computer with a JVM. 
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COINSS is written as a Java application, which means that it is executed from a 
command line prompt. COINSS and the map-based system could be converted to a Java 
applet, which would then be executed within a web browser that has a JVM, for example, 
Netscape Navigator or Internet Explorer. This would be the most convenient method of 
distributing a dynamic view of the search. This would allow searchers, family members, 
and others to have the same view of the search that the Commander has. 

COINSS uses the stochastic representation for target movement, which is 
explained in complete detail in Chapter III. COINSS and a map-based planning system 
are available by contacting Professor Gordon Bradley at NPS, Monterey, CA 93943 or 
email to bradley@nps.navy.mil. 

E. THESIS STRUCTURE 

Chapter II explains the search and rescue process and discusses how the On Scene 
Commander would use a search movement model in a search and rescue operation. 

Chapter III explains the logical search and movement elements of the Coordinated 
Inland area SAR System (COINSS). 

Chapter IV explains the graphical implementation of COINSS. 

Chapter V provides the conclusions and recommendations. 

Appendix A provides the Markov transition matrices for the five classes of targets 
defined in COINSS. 
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Appendix B provides a walkthrough of COINSS operating with a map-based 
planning system. 

Appendix C provides the code for SARgrid, and SARhex. This is the complete 
code for COINSS. 
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II. THE SEARCH AND RESCUE PROCESS 



A. INITIAL ACTION 

1. Obtaining an Initial Position 

The initial area is the last known geographical location of the target. The initial 
area can be multiple disjoint or overlapping locations, in which case a probability is 
assigned to each portion sequentially, corresponding to the probability that the target is 
located in that portion. One of three general situations for the initial area information 
usually exists [ 18 ]. 

The first is a known position. This position can come from a direct witness, or 
may even come from the target individual. The single position is used as the initial area. 

The second is a known track. In this instance the target’s route of travel is known. 
What is unknown is the position along that route of travel. The complete length of the 
known track is used as the initial area. 

The third is a known general area. This is the least favorable condition, but also 
the most likely. This instance occurs when general information regarding the habits of 
the target is known, such as where he typically hunts, or the part of the lake where she 
typically fishes. Natural boundaries, such as cliffs or rivers, can help to constrain this 
type of initial area. 

In each case, the initial areas are assigned a probability, which is the probability 
that the target is located within that particular geographic area. The combined set of 
locations creates the initial area for the search. 
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2. Defining the Search Region 



The total search region is the geographic area determined to be the most likely to 
contain the target. The bounds of this region need to account for the error in the initial 
position, the movement capabilities of the target, the amount of time the target has had to 
move, and the physical geographic constraints of the area [ 18 ]. 

The search region must be large enough to ensure that the target is located within 
it. Conversely, the total search region area should be small to minimize the probability of 
searching areas where the target is not located. Ideally, the total search region should be 
dynamically adjusted, to account for the both the target's movement, and the search 
efforts conducted. 

If available, information about the elevation and terrain within the search region 
can be used to help predict the movement of the target. 

B. SEARCH ACTION 

The search action is a looping process. The process evaluates the current 
probability information, identifies the search areas, assigns units to conduct the search, 
executes the search, and then updates the probability information based on the new search 
information. This process is continually repeated, often for a variety of search units, until 
either the target is located or the search is suspended. Because this thesis does not focus 
on the tactical decisions required to create search plans, this section will only discuss the 
role that the outcome of those decisions has on the search movement model. 
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1. Identifying the Search Area 



A search area is a sub-section of the search region that is searched by a particular 
search unit. Search areas are chosen to make best use of the abilities of the search unit in 
order to maximize the probability of detecting the target. Generally, the areas with a 
higher probability of containing the target are searched before the lower probability areas. 

2. Search Execution and Results 

After being assigned a search area, the unit searches that area. Once a search unit 
has completed its search, the unit reports its findings. If the target was found, the search 
can be concluded. If the target was not found, that information is used to update the 
probability map by Bayes Theorem. The probability that the target is in the searched area 
is reduced, and the probability that the target is in any non-searched area is increased. 

The change in probability on the probability map is directly related to the search unit's 
probability of detection of the target, given that the target is located within the area 
searched. 

Upon completion of the search result update, a decision is made whether to 
continue the search or suspend the search. This decision is an easy one when the target 
has been located. When the target has not been located, the decision takes into account 
many factors. These factors include the length of time the search has already expended, 
the probability the target will be located given that it has not been found up to this time, 
and the value of the target. If it is decided to continue, the search cycle is repeated. 
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III. LOGICAL MODEL 



A. PROBABILITY MAP COMPONENTS 

The first layer in the logical model is the probability map. This map presents the 
probability that the target is located at a certain geographic position. This probability 
follows the laws of any standard probability distribution, and computations (such as 
Bayesian updates and using Markov chains) can be applied to it. It is this property which 
allows us to update the map for both time based movement and search result information 
[ 19 ]. 

1. Elements of the Grid Square 

The logical map can be imagined as a large grid of squares. Each square is an 
individual piece that contains information relating to the complete map. Though the 
information values assigned to each square will vary from square to square, the type of 
information provided will be the same. As such, every square can be represented as an 
individual object, each having its own set of variables. The essential variables for each 
square are its position, elevation, and a velocity array. The first two of these variables are 
static. These values are set when the model is initialized, and will not change during the 
course of the search. 

The velocity matrix, V, is a four-dimensional array, with the first two dimensions 
representing the location on the probability map, and the second two dimensions being 
speed in the east-west direction and speed in the north-south direction as it would appear 
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on a common geographic map. The units for speed are feet per second. The values of the 
indices indicate the velocity in that direction. For example, V x>y ,4,3 represents the 
probability that the target is located in square S x , y moving with a velocity of 4 feet per 
second toward the east and 3 feet per second toward the north. Likewise, V x>y ,- 4,-3 
represents the probability that the target is located in square S x>y moving with a velocity 
of 4 feet per second toward the west and 3 feet per second toward the south. 

B. INITIALIZING THE PROBABILITY MAP 

There are two steps in initializing the logical probability map. The first step is to 
define the entire potential search region. The second step is to identify the potential 
initial areas for the target, as well as the respective probabilities for each area. 

Defining the potential search region is a fairly simple process. Because the 
program uses a grid array, the shape of the region will always be a rectangle. In order to 
define the bounds of this rectangle, the user provides the latitude and longitude of the 
upper left comer of the region and the lower right comer of the region. The map-based 
planning system provides these coordinates to COINSS, along with the length and width 
of the region in terms of graphical pixels. COINSS then determines the search grid array 
size. Currently in COINSS, each grid square is 16 pixels by 16 pixels on the graphical 
display. The search grid array size is then computed by dividing the search region into 
the required number of 16x16 pixel squares to cover the region. 

Defining the initial area for the target is a little more involved. A vector of 
geographic positions, generally denoted by longitude and latitude, is required. These 
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positions will define the vertices for a polygon. This polygon denotes the initial area. 

This polygon does not need to be a regular polygon, and can be non-convex. 

Each polygon is assigned a numerical percentage value, which is the probability 
that the target started in that area. The model takes this probability and distributes it 
equally among the grid squares located within the area. A grid square is considered to be 
located within the area if the center of the square is located within the area, otherwise it is 
considered to be located outside the area. 

Multiple initial areas may exist, and these areas may overlap. Each initial area is 
created in the same manner. It is conceivable that the sum of the probabilities for all the 
initial areas would be less than one. This would occur if there exists a probability that the 
target could have initially started outside of the search region. 

C. UPDATING THE PROBABILITY MAP 

Two events require the update of the probability map. The first is the passage of 
time, which requires an update due to the movement of the target. The second is a report 
from a search team, which requires an update based upon the detection/non-detection of 
the target. 

1. Time-Based Update 

Time based updates adjust the probability distribution on the map due to the 
movement of the target. These updates can occur over any time interval. Smaller time 
intervals, such as one second, come closer to simulating a continuous time movement 
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model. However, larger time steps reduce the computations, and over a large region with 
several hundred movement vectors, the loss of resolution is not significant. 

COINSS uses a constant time step, currently set to fifteen minutes. COENSS only 
computes the time-based update when prompted by the user. This allows the user to 
control the pace of the model. The time-based update consists of three major steps: 

1 . Moving the velocity array probability elements from their current square to a 
new square, based upon their velocity and the time step increment. 

2. Adjusting for changes in velocity after each probability element movement. 

3. Correcting the velocity array probability elements for elevation. 

The first step of the time-based update is to move the velocity array probability 
elements. This is done by sequentially accessing the velocity array of each square, V x , y ,jj. 
Each element of the velocity array is moved the direction and distance that is dictated by 
the velocity vector and the corresponding time step increment, Ti, to a new location, 

Sxnew.ynew- 



xnew = x + i*T } 
ynew = y + j * 7j 

The probability at this new location is evenly distributed over a square region, 
centered on the new geographical position (Figure 3-1). The probability is then assigned 
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to each grid square in proportion to the area of the probability square that is located 
within the grid square [6]. 




Figure 3-1. Movement and distribution of velocity array 
probability element 

The next step is to adjust the current velocity of the probability element. The 
velocity of the target is modeled to change according to a Markov chain. To model the 
velocity changes as a Markov chain, each velocity element is recognized as a state. The 
probability of changing velocity from state i to state inew, which can be read as changing 
from velocity i to velocity inew, is represented by a Markov transition matrix, MMj,i new 
[13]. The transition probabilities of this matrix are based on the target type of the search. 
Appendix A discusses the five target types defined in COINSS and provides the Markov 
transition matrix for each type. Figure 3-2 is the transition matrix for target type 1 , which 
is a child age 1-6. 
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Transition matrix (feet per second velocity) 
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Figure 3-2. Transition Matrix for a child age 1-6. 
Each row of this matrix sums to 1 . 



This matrix is applied individually to both the x and y components of the velocity. 
For example, from Figure 3-2, the probability of transitioning from velocity -2 (ft/sec) to 
-1 (ft/sec) is .6, and is written as MM . 2 .- 1 - The probability element, labeled V x>y> y, can 
then be distributed to the velocity array of its new square, V xneWjynew ,inew,jnew, according to 
the following: 



V. 



xnew,ynew,inew y jnew 



V, ... *MM.. * MM . . 

i r x,y y uj jyjnew 



x=xnew-T x *i, 
y=ynew—T\*j, 



For example, suppose V x>y> _ 2,-2 =.25. Then suppose that by following the example 
of Figure 1 , it is computed that 40% of the moved square lies within S xnew ,ynew, or Vyyjy 
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= .1, where x’=xnew-Ti*i and y’=ynew-Ti*j. By using the transition matrix for target 
type 1 from Figure 3-2, we obtain the resulting distribution shown in Figure 3-3. 
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Figure 3-3. Resulting distribution of V xne w,ynew,inew,jnew when 40% of 
V x,y,-2,-2 =-25 is moved to square S xne w,ynew, for target type 1 , 
small child age 1-6. The sum of all elements in this matrix is .1. 



Once the probabilities have been moved and updated for all the squares in the 
search region, the final step of the time-based update adjusts the velocity array 
probability elements to account for elevation. For COENSS, elevation can be loaded from 
a CD-ROM that contains Digital Terrain Elevation Data Level 1 (DTED1) from the 
National Imagery and Mapping Agency. This data provides the elevation in meters and 
is available to the U.S. Government only. The coverage for this data is one measurement 
every three arc seconds, or approximately every 300 feet. If elevation data is not 
available, then the elevation is considered to be 0 meters for all locations. 

The adjustments for elevation are done by again cycling through all the squares of 
the search region. At each square, the total probability of that square is computed. 

Pinit = IX„ V 

'•j 
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The next step is to cycle through the velocity array V of that square. At each 
probability element, V x . y ,j,j, of the array, a slope is computed. This slope is computed in 
the direction that the element would move during the next time step (i.e., the slope from 
S x .y to S X new,ynew , where xnew = x +i*Tj and ynew = y +j*T ( ). If the slope is too steep for 
the particular target type to easily navigate in that direction, then the probability element, 
V x.y.i.j, is reduced. If the probability is reduced, it is either reduced to 50% of its previous 
value or to zero, again depending on the severity of the slope and the target type (Figure 
3-4.). These values were determined from judgements of human tendencies. 
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Figure 3-4. Slopes at which the probability of moving is 
reduced either 50% or 100% 



After all elements of the array V in Square S x , y have been updated, the total 
probability of the square must be readjusted back to its original value, Pinit x>y , by 
computing the new total square probability, and re-scaling each element of the velocity 
array. Let V' x ,y,jj be the velocity array elements after adjustment for elevation, and let 

Pnew x, y = ZK, jj 
ij 

then 



Pin it * V' . . 



■x.y 
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2. Search-Based Update 



Each search conducted in the search region is applied to the model independent of 
any other search. The variables required for a search are the search unit name, the search 
area, and the probability of detection. The search unit name is a label to help identify this 
particular search instance. The search area is a vector of geographic positions, generally 
denoted by longitude and latitude. These positions will define the vertices for a polygon. 
Unlike the search region, the search area is not required to be a rectangle and it may be a 
non-convex polygon. The probability of detection is a value between 0 and 1 that is 
provided by the user. 

Search based updates adjust the probability distribution for all the squares in the 
logical map due to search unit results. These updates use Bayes theorem to account for 
unsuccessful search results. This result is computed by using conditional probability rules 
and Bayes theorem. The following two events are defined: 

D xy : The target is detected at location S x , y by the search effort allocated to S x , y 

T x>y : The target's current location is S x , y 

Then the probability of detection for the search unit, conditioned on the event that 
the target is located in the searched area, is defined as: 

P[D x , y I T x , y ] 

and the probability of non-detection, conditioned on the event that the target is located in 
the searched area is defined as: 

P[D’x, y I T x , y ] = 1 - P[D x , y I T x , y ] 
where the superscript denotes negation. Likewise: 

P[D xy \T xW ] = Oy(x,y)*(x',y') 
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S x , y is determined to be located within the search area if and only if the center of 
square S x<y is located within the search area polygon. 

When a search party reports that they have searched an area and did not find the 
target, the model is updated. Each square in the search region must first be updated to 
show the probability that the target is located in that region and was not detected by the 
search party. This is defined as: 






Using Bayes theorem, we arrive at a posterior probability that is the probability of 
the target being located at location S X ',y, given that it was not detected at location S x>y . 
This is defined as: 



P[T,AP[D'\T,,] 

PIT . . I D' 1 = ° x -> — 1^-i. Vx v x' v' 

1 



It is worth noting that all squares in the search region are being updated, not just 
the squares located within the search area. For all the locations S x -y located outside the 

search area, P[D' X y |T X . y . ] = 1 , and therefore P[T x . y . f| D x y ] = P[T X . y . ] . These values are 

used in the above equation to create a new probability distribution for the location of the 
target. 
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IV. VISUAL REPRESENTATION OF THE MODEL 



A. COMPONENT-BASED DESIGN 

COINSS was designed to be a component-based system. COINSS provides the 
algorithms, methods, and logic for executing a S AR mission plan. It is designed to take 
advantage of well-defined reusable components from existing map-based planning 
systems, such as the ability to read latitude and longitude or to provide overlay displays. 

The "Loosely Coupled Components" project designed by Dr. G.H. Bradley and 
Dr. A.H. Buss, both of the Naval Postgraduate School, provides one such map-based 
planning system [4]. Their project creates planning systems by integrating various 
components as required. Components can be maps, algorithms, data overlays, or any 
other item necessary to complete a task. The nature of the system is such that the planner 
only needs to assemble the distinct components needed at that time. If a need for a new 
component arises during a planning operation, the planner can obtain or implement that 
component and add it to the planning system without having to restart the entire 
operation. It also allows designers to concentrate on creating algorithms or methods to 
solve problems, and not continually "re-invent the wheel" by developing a graphical 
interface for their specific problem. 

B. BASIC CONSTRUCTION OF THE VIEW 

1. Grid Dimensions 

To present a COINSS probability map on a visual display, a translation from the 
squares on the logical probability map to the pixels on the monitor is required. A 
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standard method of presenting probability maps on a visual display is by using different 
colors or shades of a single color to represent the probabilities. The visual display square 
size needs to be large enough to allow visual discrimination of color and the squares must 
be small enough to provide useful search information (i.e., a search region with only two 
squares would not provide much meaningful search information). Finally, the logical 
array needs to be small enough to provide movement and search updates in a reasonable 
amount of computation time. COINSS was developed on a PC with a 166 MHZ CPU. A 
16x16 pixel square is a reasonable choice for this computer. 

2. Probability Display 

COINSS uses 10 shades of the color red to discriminate among the different 
quantitative probability values. The red shades are equally spaced, which is achieved by 
manipulating the RGB color value of the pixels. The RGB value ranges from 0 to 255. 
The shades were created by modifying the values in increments of 25 from 5 to 255. Red 
was chosen as the color to use, and the darkest red was achieved with RGB values 
2551010, while the lightest shade was achieved with RGB values 25512251225 [8]. 

The variable maxSquareProb contains the value of the largest probability in the 
search region. 

maxSquareProb = max{^ V^.jh^y 

Each square is assigned a color, which is based on the probability of the square 
compared to the maxSquareProb. The maxSquareProb is assigned the purest color value 
(2551010). Squares containing a probability that is 90%-99% of maxSquareProb receive 
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the next highest purest value (255125125). This continues at 10% intervals down to the 
squares that have probabilities that are between 10%-19% of the maxSquareProb. These 
squares receive the least pure color value (25512251225). The squares with probability 
less than 10% of the maxSquareProb are assigned to be completely transparent, and 
appear as blank squares on the probability map. 

COINSS also makes use of opaqueness. The opaqueness of a cell determines how 
well you can see through the cell on the top layer and view images on the layer(s) below 
it. An opaque value of 255 is completely opaque and does not allow the viewer to see 
any part of the layer below it, while an opaque value of 0 is completely transparent [8]. 
COINSS currently assigns each colored square in the probability map an opaqueness 
value of 190. This allows the layer immediately below the probability map to be seen, 
with little loss in color discrimination. 

C. INTERACTING WITH THE VIEW 

1. Selecting a Geographic Map 

COINSS is designed to operate in conjunction with a map-based planning system, 
such as the Loosely Coupled Components system created by Bradley and Buss [4], The 
map-based system is responsible for obtaining a geographic map from any available 
source. The planning system is also responsible for providing an interface to COINSS, 
specifically, a reference to the latitude and longitude of the locations on the map, as well 
as a representation of the scale of the map provided in pixels per unit longitude. The map 
can have any geographic or topological features or symbols. This map is merely a visual 
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reference for the planner, and will have no impact to the logical computations that take 
place in the model. 

2. Defining the Search Region 

Identifying the bounds for the search region is the next step. This is a user- 
defined region that is provided by the map-based planning system. The search region is 
required to be a rectangular region, and is defined by two points, the upper left comer and 
the lower right comer. A simple method is to allow a mouse click at the location for the 
upper left comer, followed by another click at the location for the lower right comer. 
These two values are then provided to COINSS to define the search region. 

3. Defining the Initial Areas 

Initial areas are also user-defined areas that are provided by the map-based 
planning system. Again the mouse is used to outline the initial areas. As discussed 
earlier, a vector of geographic positions is required. These positions define the vertices 
for a polygon. The mouse clicks define each vertex, in order. Once the polygon is 
completed, it is assigned a numerical value, which is the probability that the target started 
in that area. COINSS takes this probability and distributes it equally among the grid 
squares located within the region, and across all the possible velocity vectors. 

This is an iterative process, conducted on each initial area in a sequential manner. 
Each subsequent initial area probability is added to the existing probability map. If two 
or more initial areas overlap, then the probability associated with the overlapped area is 
equal to the sum of the individual probabilities for that square. 
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4. Defining the Search Area 



Search areas are the third case of user-defined areas. These are also provided by 
the map-based system, and are created in the same manner as the initial areas. The 
numerical value assigned to each search area is the probability that the unit searching that 
area will locate the target if the target is in that area. This probability is used to apply 
Baysian updates to the probability map. 
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V. CONCLUSIONS AND RECOMMENDATIONS 



COINSS provides an inland search and rescue model which implements target 
movement during the search execution. Target movement is modeled as a discrete time 
Markov chain, based upon the type of the target. Geographic elevation is also accounted 
for in the movement model by limiting movement in areas with large changes in 
elevation. 

COINSS is a component-based system, designed for implementation with a map- 
based planning system. This design demonstrates that a component-based design can 
reuse components. Allowing different component systems to use the same map-based 
system interface increases the efficiency of the hardware, reduces user training, and 
reduces software cost and redundancy. 

This thesis provides several areas for future research. The movement model is 
very rudimentary. The COINSS Markov movement model is based upon historical data 
that only provides the total distance traveled from initial datum [9] [17]. More detailed 
study of target movement rates and tendencies would improve the movement model. 

Integrating additional information regarding the terrain would also improve the 
movement model. Identification of impediments such as heavy brush, rivers, and lakes, 
or easier traveling terrain such as trails, paths, and roads would improve the realistic 
modeling of movement, and create a more realistic probability model. 
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APPENDIX A. MARKOV MOVEMENT MATRICES 



At each movement phase update, the velocity vector of the target changes as 
discussed in Chapter III. To implement the movement model using a Markov chain, a 
discretized version of a Kalman filter was created [22]. 

Using a discrete version of the Omstein-Uhlenbeck (OU) process [22], velocity is 
modeled as 

X t = cX , + W , 

where X , is the velocity of the target at time t, W t is an independent identically 
distributed normal random variable with mean 0 and variance q, and c is a constant. The 
values for c and q determine the variance of the velocity for the target. Furthermore, 

A 

c = e T 

where A is the time interval at which updates are made, and 1/x is the expected rate (in 
time units) that the target changes velocity [22], A is an artificial variable which is 
determined by COINSS, while t is a characteristic of the target that is estimated from the 
data. 

The variance of the total movement of the target can be expressed as [19] 

1 _£ 

Var(X(5)) = -(2s 2 T)[25-r[3-4e ' +e r ]] 
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2 

where s' is a target characteristic which identifies the variance of velocity, and 5 is the 



time interval over which the target has moved. By computing the variance for different 



target types from Syrotuck [17] and Hill [9] (Figure A-l.) and estimating values for 8 and 
x, s' can be solved for, and then q can be computed from [22]: 

2/i 2\ 

q = s (i-c ) 
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Figure A-l. Variance of total movement (in miles) for found targets 
in Hill[9] and Syrotuck[17], 



The data provided by Hill and Syrotuck in Figure A-l have two shortcomings. 
First, the data is only for found targets. As such, we must assume that all targets move in 
a fashion similar to that of the found targets. Second, the data has no time associated 
with it. The longer a target remains undetected, the larger the variance of his movement. 
However, this is the best data available on lost targets. Therefore, the data was used, 
along with a generic view of human movement patterns to create one model of 
movement. 
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After c and q have been estimated, the Markov transition matrix can be 
constructed using the assumption that velocity changes will vary according to a standard 
normal distribution. Figure A-2 provides a chart of all the values for x, which were 
estimated based upon human tendencies, and s 2 , which was computed with an estimated 5 
of 24 hours. Figures A-3 through A-7 are the transition matrices for each target type 

when A = 15 minutes. The indices for the transition matrices are in feet per second. The 

© 

target types are [9]: 

1 : children under age 6, 

2: children age 7-12, 

3: despondents/walkaways, 

4: hunters, 

5: hikers/adults. 

Target Type x(sec) S 2 (ft 2 /sec 2 ) 

1 
2 

3 

4 

5 

Figure A-2. Computed values for target movement 
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Figure A-3. Transition matrix for target type 1 
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Figure A-4. Transition matrix for target type 2 
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Figure A-5. Transition matrix for target type 3 
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Figure A-6. Transition matrix for target type 4 
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Figure A-7. Transition matrix for target type 5 
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APPENDIX B. COINSS SCENARIO 



This is a walk through of COINSS implemented on a map-based planning 
system. The system used to implement COINSS in this scenario is Special 
Operations Forces/Loosely Coupled Components (SOFLCC) [2], 
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Figure B-l. The standard map-based planning system. The latitude and longitude of the 
mouse position is displayed in the lower right of the display. 
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Figure B-2. To add a new overlay to the map, click the “Add a new Overlay” Icon. This 
then provides the option to add a (new) blank layer, or add a (pre-existing) layer from a 
file. 
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Figure B-3. The search region is created by clicking once on the map at the north-west 
comer of the search region, then double-clicking at the south-east comer of the region. 
The map has been faded to provide a clearer view of the search region. 
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Figure B-4. Select the Show Overlay icon to select the search region for initialization. 
Select the Run Algorithm option to execute a SAR algorithm. 
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Figure B-5. Select the Add Algorithm option to dynamically load the SAR algorithms. 
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Figure B-6. Select the SAR.initialRegion to create the search region. 
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Figure B-7. The initialRegion algorithm requires a name as an input parameter. 
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Figure B-8. The system queries for the type of search target, which is contained in a file. 



46 












class m l, a r Ttf.tr s c Konl q. Graph ££ 
Map laiy'er Lpoer Uti* tv35' 16’C£ 



oublic stetic vo d aofcc. algorithm G^R.initiaRecicn^milarmy.trac.l-ioiia.Ciraohjava.Iang.Ct'ing) 
public st?tic vo d eof ct. algorithm SAR.initialArea(m l.anij<.1r5C.kgnig.Graphjav/a.l3ng.Btrjng,javalang.S1rin£) \ 
oublic Etstic vo d sot cc. algorithm y^K.S0archAr3a(rnil.armv.lrac.koniq.Cjraprija ( /s.larq atnnq.iava IsnuiJlrinq) \ 
oublic static vo d sofct. algorithm SAR.mo^ementirnil amytrac.Horiig.Oiapl-.java lang.9trirg java lanerEtirg) 
oublic Etstic s*o d eofcc. algorithm S4R.sleepUpdate(ni .srmv tra:.konig.Graph,java.lang.S1rinc,java.lang.Strirg«5 
nuhlir Fistic vn rt <nf rt algorithm RAR rp^riFtRvatihnD#a(nt Srtny tta* kmiq Graph, java tanrj R1rin$ I 

oublic static vo d sofcc. algorithm SAR.updaleProblmaceO 



Show L«j 



Hide Layer 









Add Algorithm 



Run Algor Khm 



Close 



l^tiowuwrtws 



uwyimnim 



» £A *i w! 









Dvu«i 

CJW082 

LJV083 

; * ivos* 

LJV085 

illwoo: 



•i 

^jw5»ib 
i l - jwo^i 

t ?W0 7 2 
JW3 
„ JW0 7 4 

s.y >- 

rtegMfte. [*** 

Flw |il Flet | K ! ) 



JWUft 

_VQ7$ 

^.iW077 

W078 

IJWCF9 

liwooo 



□V368 

;.JW3?9 

::: iwjso 



v* 



V^SWftYrtrwlwiVAVtSV.' 



r >m 






jN3$* 10*026,*“; pow *055\i 



Choo;c DTCO file 



Figure B-9. If elevation data is available, the user can import that data into the system. 
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Figure B-10. Add a new overlay to create the initial area. The area is created by single 
clicks of the mouse at each location of the map which is a vertex for the area. Double 
click at the last vertex to complete the area. 
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Figure B-l 1 . Use the Show Overlay icon to select the initial area graph, and run the 
initialArea algorithm. This algorithm requires a name and a probability between 0 and 1 
as parameters. 
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Figure B-12. The probability map updates to show the current probability distribution. 



50 











$mtcom%mf^n^ 

class roil.ainTtrat.tomg.Oraph@: 
Map layar Mpp*r Lat~ N35* 16 1 C* 



Si 



Atftf Alyaritliin 



|0| Show Overfcys 



i: 



W-t 



Run AlQorrttun 






M*S*>' 



i >4^t ! ’ , ' ; i 



• A v\ *• "S 



mviw&w-iwtm'wm# 



:v << 4fvs&* 



Running Algorithm 



c&*s. miimm:i^&:. & %%H : : 

public ctati: void colic: algcrrhTt GAR JnltiolRDgion(mll.arnv.ttac.konig G ophjavo long String) 
pnhlic stath vnlrl antin' algrrrhTi FW.R initwIAma^il 3rrny ’rar knrlc Oratrh jFva lang ftr n^ava lang S*rhg) 
public statl: void so(lc:.algcrl:lmSAR.secn:h. i \rea(m l.am/.lrcC.konlg.Oraphjava.lang.strlngJavalang.sirlnc) 
public statb void soitc:;algcrthri.GARjroveire^rni!.arrnytr9c koilg.CJia^Java.Iang.OHngJavaJanc.Otrfng) 
public statl: void cotlc:.algcrl-h"n.SAR.'Slee[;IJpd?t6(mjl.ariny.t , ac.<cn g Gr?p^avalang.S:riTgj3\>a.lang. Sting) 
public stati: void sonc:.aiqcrr.nTibAK.rea3biev3io'iusK<niii.aririvt'ac.<cn a Grspi, ava.iana.y:rnq) 
public stati: void soflc:.algcrl:hii.SAR.upca:e :> toblniageO 



Figure B-13. Use the Show Overlay icon to select the search region, and run the 
movement algorithm. The algorithm requires a name and a length of time in minutes. 
The probability map is updated based on the target movement that would occur during 
that time interval. A second movement option is available, the sleep update. This option 
also requires a name and a length of time in minutes. The time on the probability map is 
updated, but no target movement occurs, because the target is assumed to be sleeping. 
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Figure B-14. Add a new overlay to create the search area. The area is created by single 
clicks of the mouse at each location of the map which is a vertex for the area. Double 
click at the last vertex to complete the area. 
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Figure B-l 5. Use the Show Overlay icon to select the search area graph, and run the 
searchArea algorithm. This algorithm requires a name and a probability between 0 and 1 
as parameters. The probability map is updated based on the search results. 
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Figure B-16. Double clicking on the probability map displays the current date and time of 
the search. 
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APPENDIX C. COINSS CODE 



This appendix contains the code for SARhex, and SARgrid. In addition, it 
discusses the driver program necessary to interface COINSS with a map-based planning 
system. SARhex is the basic building block for COINSS. SARgrid contains all the 
methods which operate the model. 

The COENSS software is connected to a map-based planning system through a 
driver program that invokes the methods in COINSS. The driver invokes these methods 
by obtaining the necessary information from the map-based planning system. The driver 
then converts the information into the correct parameter form to call the methods in 
COINSS. The driver then invokes the methods using the proper parameters, and provides 
the results back to the map-based system. 

To initialize the search region, the driver instantiates an object “SearchRegion” of 
class SARgrid. The constructor requires the latitude and longitude of the northwest 
comer and the southeast comer, each being double variables with units degree. The 
constructor also requires the height and width of this region on the screen in pixels, each 
being int variables. 

To specify the initial area(s) the method problnitialDatum is invoked on the object 
SearchRegion. The first parameter is a two-dimensional array of type int. Each column 
of the two-row array is a vertex of the polygon that defines one initial area. Each column 
is the location, in grid squares (x,y), of a point. The second parameter is a double 
variable that is the probability (0.0<j < 1.0) that the target started in this area. 
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To specify the search area(s), the method problnfoUpdate is invoked on the object 
SearchRegion. The parameters are the same as the method problnitialDatum, where the 
double variable is the probability the search team will find the target, given he is located 
in the given search area. 

To specify movement, the method probTimeUpdate is invoked on the object 
SearchRegion. This method requires no parameters. It will update the probability map 
for movement of one time step, based upon the constant time step defined by S ARgrid 
(currently 15 minutes). 

To specify time passing during the night, the method updateTime is invoked on 
the object SearchRegion. The only parameter is a double variable. This double is the 
number of seconds that have elapsed while the target sleeps. 

To present the probability map, the method getProbMapIcon is invoked on the 
object SearchRegion. This method requires no parameters. It returns the current 
probability map as an Icon variable. 
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public void resetTempProb ( ) { 



public void incTempProb ( int i, int j , double value) { 
tempProb [ i ] [ j ] + = value ; 
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public double getProb(int 
return (prob [ i ] [ j ] ) ; 
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public static int Scale = 16; //Dimension of each grid hex in pixels 



private int targetType; //Numeric value for the type of 

//target being searched for 
//I - small children (1-6) 

112 - children (7-12) 

// 3 - despondents/walkaways 
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private short [][] cellDTEDl; //Array of elevation data for a 1 degree 



when created, SARgrid requires the latitude and longitude of the northwest corner 
(latl,longl) and southeast corner ( lat2 , long2 ) , as well as the number of pixels in 
the region, both vertically (pixl) and horizontally (pix2) 
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* * * * 



cellDTEDl = null; 



// instance methods 
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for(int x=0;x<=gridWidth-l;x++) 

for(int y=0 ; y<=gridHeight-l ;y++ ) { 

if (inPoly (ini tialAreaAr ray, x,y) ) { 
square [x] [y] . setProb (prob/count ) ; 



// total the probability in the entire search region 
// then normalize each square's probability, so that the 
// total search region probability = 1 
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if (newGridI<=0) {incOuter (3 , square[i] [ j ] . getProb (x, y) ) ; 
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public void problnf oUpdate ( int [ ] [ ] searchAreaArray , double pDetect) { 



double totalProbStart ; //total prob of target in search region 

double totalProbEnd; //prior to update, and after update 
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Modified code from Inpoly, C 

Copyright (c) 1995-1996 Galacticomm, Inc. Freeware source code. 
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xt) == (xt <= xold) /* edge "open" at left end */ 
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for ( int j =0 ; j <=gridHeight-l ; j ++ ) 

for(int i=0; i<=gridWidth-l; i++) { 
value = computeElevColor ( i , j ) ; 
for (int x=Scale* j ; x<=Scale* j +Scale-1 ; x+ + ) 
for ( int y=Scale*i ;y<=Scale*i+Scale-l ;y++) 



printArray [x*Scale*gridWidth+y ) = value; 
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new double [maxVel ] [maxVel] ; 



for (int i=0 ; i<=maxVel~l ; i++) { 
for (int j=0; j <=maxVel-l ; j++) { 

moveMatrix [ i ] [ j ] =in . readDouble ( ) ; 
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f ileDialog . setMode ( FileDialog . LOAD) ; 
f ileDialog . show( ) ; 
fileName = f ileDialog . getFile () ; 
if (fileName == null) { 

fileName = "dted/w098/n3 1 . dtl " ; 
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System, out . print (gridHeight+ "x" +gridWidth) ; 



horPixels = gridWidth*Scale ; //the scaled search map encompasses , 
lrLat = ulLat-incLat*gridHeight ; //assigns a new lower right Lat/Long 
lrLong = ulLong+incLong*gridWidth; //based on this new scaled coverage i 
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This getter returns the elevation at a given longitude/ latitude 



The method assumes that the long/lat provided as (i,j) are within 
the same degree of long/lat selected during the readDTEDl method 









O 




















O 




















CM 




















T 1 




















* 








































x 




















C 




















•r| 




















■ 




















1 

*n 






























O) 










— ' 










G 










, — . 










O 










X 










X 










G 










O 




















G 










' — 










X 










1 — 1 










\ 




















tn 




















G 










O 










O 










O 










X 




•O 






CM 


T3 








X 










X 


G 








3 




0 






* 


0 












t — 1 






. . 


U 








« — 1 




X 






X 


(U 








03 




3 






, — . 


CO 












O 






X 




















G 


X 


















-H 


G 








X 




«. 






— 


•H 


X 






G 




*H 






1 




03 






3 












** 


> 






O 




0 






— 


X 












1 — 1 






— 


CO — 


0 










X 






1 


u — 


1 — 1 






X 




3 






O 


X T3 


X 






X 




O 






O 


x C 


3 






03 




*0 






CM 


O 


O 




»«. 


2 




— " 






X 


X 0 






X 


' — 




G 






— - 


G 0 


— - 




X 


„ — , 




O 






- — . 


*H CO 


0 




T3 


X 




X 






X 


— ' > — ■ 


3 




X 


G 




X 






G 


a) — 


X 




— £ 


X 




03 








u x 


03 




— T3 


' — - 




> 








03 CO 


> 




O) x 






0 






* — < 


3 u 


1 — 1 




G U 


II 




1 — 1 






t— 1 


D 1 X 


0 




O tn 






W 






Q 


10 X 


X 




X 


X 


... 


x 






W 


X “ 


-H 


. . ^ 


U II 


0 


. 


0 






&H 


0 0 


a 


tn 0 


X 


X 


1 — 1 


Cn 






Q 


Ol U 


X 


G 


A X 


X 


0 




•«. 




H 


03 


x — 


O 11 


X 0 


^ a 


X 


0 


— .. 




rH 


X 3 


0 1 — 1 


X 


03 X 


0 X 


X 


1 — 1 


X 




0 


CD O 1 


tn 0 


1 — 1 1 — 1 


> X 


CO 


a 


X 


- — 


- — 


0 


X CO 


X 


3 0 


— a 


1 — 1 


X 


3 


CO 


CO 


— 


Dh — 


X X 


V X 


X X r- 


- 0 r~> 


— * 


O 


X 


X 


G 


< G 


G a 


X X 


X 




G 


*0 


03 


03 


U 


W Vh 


x x 


03 a 






U 




• 


• 


3 


3 




j> X 


0 




3 


U 


X 


X 


X 


U X 


0 X 




CO 




X 


- H 


X 


X 


0 


*H 0 


X c 


X 


1 — 1 




0 


1 — 1 


03 


03 


U 


X H 


1 — 1 *H 


X 


0 




u 


x 


£ 


2 




X 


X 










3 


II 


II 




3 


3 










a 


•H 


*n 


r~~% 


a ^ 


a 











* * * 



public int getYPixelValue (double val){ 
int yPixel; 



yPixel = (int) (Math, round ( (ulLat-val ) /incLat) ) ; 
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